PROCEEDINGS 

OF SCIENCE 



J/m at high temperatures in anisotropic lattice QCD 



in 
o 
o 

(N 

Or 
<D 

(N 



> 
(N 
OS 

o 
in 
o 

D 



X 



Hideaki lida and Noriyoshi Ishii 

Department of Physics, H-27, Tokyo Institute of Technology, 
Oh-okayama 2-12-1, Meguro, Tokyo 152-8551, Japan 



E-mail: iida@th.phys.titech.ac.jp ishii@th.phys.titech.ac.jp 



Takumi Doi 

RIKEN BNL Research Center, Brookhaven National Laboratory, 
Upton, New York 11973, USA 
E-mail: |doi@ quark . phys . bnl . gc 



JOV 



Hideo Suganuma 

Department of Physics, Kyoto University, 
Kitashirakawaoiwake, Sakyo, Kyoto 606-8502, Japan 



E-mail: suganuma@ scphys . kyoto-u .ac.jp 



J/*¥ and r\ c above the QCD critical temperature T c are studied in anisotropic quenched lattice 
QCD, considering whether the cc systems above T c are compact quasi-bound states or scattering 
states. We adopt the standard Wilson gauge action and (a) -improved Wilson quark action with 
renormalized anisotropy a s /a t = 4 at j3 = 6.10 on 16 3 x (14 — 26) lattices, which correspond to 
the spatial lattice volume V = L 3 ~ (1.55fm) 3 and temperatures T ~ (1.11 — 2.07)T C . To clarify 
whether compact charmonia survive in the deconfinement phase, we investigate spatial boundary- 
condition dependence of the energy of the cc systems above T c . In fact, for low-lying cc scattering 
states, there appears a significant energy difference AE = E'(APBC) — £(PBC) between periodic 
and anti -periodic boundary conditions as AE ~ 1yfn? c + 3n 2 /L 2 — 2m c (m c : charm quark mass) 
on the finite-volume lattice. In contrast, for compact charmonia, there is no significant energy 
difference between periodic and anti-periodic boundary conditions. As a lattice QCD result, we 
find almost no spatial boundary-condition dependence for the energy of the cc system in J 'f}f and 
tj c channels for T ~ (1.11 — 2.07) T c , which indicates that J and r\ c would survive as compact 
cc quasi-bound states below 2T C . 
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1. Introduction 

Since QCD is established, the quark-gluon-plasma (QGP) phase has been studied with much 
attention as a "new phase of matter" at high temperatures both in theoretical and experimental 
sides Ijjjimfll. In recent years, QGP creation experiments are actually performed at SPS [J3j] and 
RHIC [Q] in high-energy heavy-ion collisions. As an important signal of the QGP creation, 7/*P 
suppression [Q, |2|] was theoretically proposed and has been tested in the SPS/RHIC experiments. 
The basic assumption of J /*¥ suppression is that //*P disappears above T c due to vanishing of the 
confinement potential and the Debye screening effect [Q]. 

Very recently, some lattice QCD calculations indicate an interesting possibility that J ' f*¥ and 
tj c seem to survive even above T c [Q, ^, [7|], which may lead a serious modification for the 7/*P 
suppression scenario in QGP physics. However, as a possible problem, the observed cc state on 
lattices may not be a nontrivial charmonium but a trivial cc scattering state, because it is difficult 
to distinguish these two states in lattice QCD. 

In this paper, we aim to clarify whether the cc system above T c is a compact quasi-bound 
state or a scattering state, which is spatially spread. To distinguish these two states, we investigate 
spatial boundary-condition dependence of the energy of the cc system by comparing results in 
periodic and anti-periodic boundary conditions. If the cc system is a scattering state, there appears 
an energy difference AE between the two boundary conditions as AE ~ 2-y/ m 2 . + 3n 2 /L 2 — 2m c 
with the charm quark mass m c on a finite- volume lattice with L 3 . If the cc system is a compact 
quasi-bound state, the boundary-condition dependence is small even in finite volume. In Ref.[Q], 
this method is actually applied for distinction between a scattering state and a compact resonance. 

2. Method to distinguish a compact state from a scattering state 

To begin with, we briefly explain the method to distinguish a compact state from a scattering 
state in term of its spatial extension. For this purpose, we investigate the cc system in the periodic 
boundary condition (PBC) and in the anti-periodic boundary condition (APBC), respectively, and 
examine spatial boundary-condition dependence for the cc system. Here, in the PBC/ APBC case, 
we impose periodic/anti-periodic boundary condition for c and c on a finite- volume lattice. 

Table 1: Periodic boundary condition (PBC) and anti-periodic boundary condition (APBC): the relation 
between spatial boundary condition and the minimum momentum \p m i n \ of c, c and compact charmonia cc. 
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For a compact cc quasi-bound state, the wave function of each quark is spatially localized and 
insensitive to spatial boundary conditions in lattice QCD, so that the charmonium behaves as a 
compact boson and its energy in APBC is almost the same as that in PBC For a cc scattering 
state, both c and c have non-zero relative momentum p m { n = (±£,±£,±£), i.e., \p m [ a \ = V3n/L 



184/2 



J ' /m at high temperatures in anisotropic lattice QCD 



Hideaki Iida 



in APBC, while they can take zero relative momentum p m { n = in PBC. In fact, if the cc system is a 
scattering state, there appears a significant energy difference AE between PBC and APBC due to the 
finite lattice volume of L 3 , and it is estimated as AE ~ 2y m 2 + 37T 2 /L 2 — 2m c . In our lattice QCD 
calculation, \p m i n \ and AE for the cc scattering state are estimated as \p m i n \ = ~ 0.69GeV 

and AE ~ 2 v / m 2 + 3k 2 /L 2 - 2m c ~ 0.35GeV for L ~ 1.55fm and m c ~ 1.3GeV. 



3. Anisotropic lattice QCD 

In this paper, we adopt anisotropic lattice QCD for the study of high-temperature QCD. In 
lattice QCD at temperature T, (anti)periodicity is imposed in the temporal direction with the period 
l/T, and hence it is technically difficult to measure temporal correlators at high temperatures. To 
overcome this problem, we use the anisotropic lattice with anisotropy a s /a t = 4. Owing to the finer 
temporal mesh, we can obtain detailed information for temporal correlators. 

For the gauge field, we adopt the standard plaquette action on an anisotropic lattice as []|, ^] 

S G = t^r I ReTr{l-^ 7 (,)} + A 7G £ ReTr {l-/>. 4 ( 5 )}, (3.1) 

1G s ,i<j<3 iV c j,j<3 

where P^ v denotes the plaquette operator. In the simulation, we take j3 = 2N c /g 2 = 6.10 and the 
bare anisotropy y G = 3.2103, which lead to renormalized anisotropy as a s /a t = 4.0. The scale is set 
by the Sommer scale Tq 1 = 395MeV. Then, the spatial and temporal lattice spacing are estimated 
as aj l ~ 2.03GeV (i.e., a s ~ 0.097fm), and a" 1 ~ 8.12GeV (i.e., a t ~ 0.024fm), respectively. The 
adopted lattice size is 16 3 x (14 — 26), which corresponds to the spatial lattice size as L ~ 1.55fm 
and the temperature as T = (1.11 — 2.07)T C . We use 999 gauge configurations, which are picked 
up every 500 sweeps after the thermalization of 20,000 sweeps. 

For quarks, we use 0(a)-improved Wilson (clover) action on the anisotropic lattice as [||, ^] 

x,y 

K(x,y) = 8 x>y - K t {(\ - Y4 )U 4 (x)S x+ly + (1 + y A )u\{x - 4)8 x _ i y } 
~ KsY,{(r-YdUi(x)S x+ly } + (r+n)Uj(x - i)8 x _-J 

i 

- KiCEj^Oi4Fi48^y-rK s C B Y, (J ij F ij S x,y, (3.2) 

i i<j 



which is anisotropic version of the Fermilab action QIOQ. K s and k, denote the spatial and temporal 
hopping parameters, respectively, and r the Wilson parameter, ce and cb are the clover coefficients. 
The tadpole improvement is done by the replacement of Uj(x) — ► Uj(x)/u s , U 4 {x) — > U 4 {x)/u t , 
where u s and u t are the mean-field values of the spatial and the temporal link variables, respectively. 
The parameters K s , K t ,r,CE,CB are to be tuned so as to keep the Lorentz symmetry up to 0{a 2 ). At 
the tadpole-improved tree-level, this requirement leads to r = a t /a s , ce = 1/ {u s u 2 ), cb = 1/ m 3 an d 
the tuned fermionic anisotropy Jp = (u t K t )/(u s K s ) =a s /a t . For the charm quark, we take fc = 0.112 
with \/k = 1/ (u s K s ) — 2{yf + 3r — 4), which corresponds to the hopping parameter in the isotropic 
lattice. The bare quark mass mq in spatial lattice unit is expressed as m = |(i — 8). We summarize 
the lattice parameters and related quantities in Table 2. In the present lattice QCD, the masses of 
//*P and T] c are found to be mj^ ~ 3.07GeV and m T , i , ~ 2.99GeV at zero temperature. 
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Table 2: Lattice parameters and related quantities in our anisotropic lattice QCD calculation. 

/3 lattice size aj l a^ 1 Yg u s u t jp K 

TO 16 3 x (14 - 26) 2.03GeV 8.12GeV 3.2103 0.8059 0.9901 4.0 0.112 



4. Temporal correlators of cc systems at finite temperature on anisotropic lattice 

To investigate the low-lying state at high temperatures from the temporal correlator, it is prac- 
tically desired to use a "good" operator with a large ground-state overlap, due to limitation of the 
temporal lattice size. To this end, we use a spatially-extended operator of the Gaussian type as 

0(t,x) = iV£exp | c(t,x+y)rc(t,jt) (4.1) 

in the Coulomb gauge N is a normalization. The size parameter p is optimally chosen in 

terms of the ground-state overlap. T = y k {k=\—3) and r = y 5 correspond to 1 (J/^V) and (rj c ) 
channels, respectively. The energy of the low-lying state is calculated from the temporal correlator, 

G(0 = ^L<O(?,x)O + (0,0)), (4.2) 



X 



where the total momentum of the cc system is projected to be zero. 

In accordance with the temporal periodicity at finite temperature, we define the effective mass 
w*eff (0 from the correlator G(t) by the cosh-type function as [JTTt] 



G(t) _ cosh[m eS (t)(t-N t /2)} (4J) 



G(t + l) cosh[meff(f)(f + l-iVf/2)] 

with the temporal lattice size N t . In the plateau region of m e ff(f), m^{t) coiTesponds to the energy 
of the low-lying cc state. To find the optimal value of p, we calculate the correlator G(t) for p = 
0.2,0.3,0.4 and 0.5fm at each temperature, and examine the ground-state overlap by comparing 
G(t)/G(0) with the fit function of g &t (t) =Acosh[m(t — N t /2)] |[n|]. As a result, the optimal size 
seems to be p ~ 0.2fm for cc systems. Hereafter, we only show the numerical results for p = 0.2fm. 

5. Lattice QCD results for the cc system above T c 

We investigate the cc systems above T c both in J/ X ¥(J P = 1~) and 7] C (J P = 0~) channels. For 
each channel, we calculate the temporal correlator G(t) and the effective mass ra e ff(?) defined by 



Eq.(|4.3|) both in PBC and APBC, and examine their spatial boundary-condition (b.c.) dependence. 

Figures 1-4 show the effective-mass plot m e g(?) of the cc system in the J pH channel for 
p=0.2fm. From the cosh-type fit for the correlator G(t) in the plateau region of m e ff(f), we extract 
the energies, £(PBC) and ^(APBC), of the low-lying cc system in PBC and APBC, respectively. 
Table 3 summarizes the cc system in the //*P channel in PBC and APBC at each temperature. 

As a remarkable fact, almost no spatial b.c. dependence is found for the low-lying energy of 
the cc system, i.e., AE = ^(APBC) — .E'(PBC) ~ 0, which is contrast to the cc scattering case of 
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Figure 1: Effective mass of J at T = 1 . 1 1 T c in 

lattice unit. The circular/triangle symbols denote Figure 2: Effective mass of J/*¥ at 
the results in PBC/APBC. The dashed/dotted lines The symbols, lines, units are same as 
denote £(PBC/APBC) in PBC/APBC. 



T = 132T C . 
Figure 1. 




CD 
3= 
CD 




Figure 3: Effective mass of// 1 ? at T = 1.617;. Figure 4: Effective mass of J at T = 2.Q1T C . 
The symbols, lines, units are same as Figure 1 . The symbols, lines, units are same as Figure 1 . 



AE ~ lyjrn^. + 3n 2 /L 2 — 2m c ~ 0.35GeV for L ~1.55fm and m c ~ 1.3GeV as was discussed in 
Sect.2. This result indicates that J ' [*¥ survives for T = (1.11 — 2.07)T C . 

Table 4 summarizes the cc system in the T] c channel in PBC and APBC at each temperature. 
Again, almost no spatial b.c. dependence is found as AE = /i(APBC) — £'(PBC) ~ 0, and this 
result indicates that 7] c also survives for T = (1.11 — 2.07)T C as well as Jf¥. 

In contrast to J /*¥ and r\ c , our preliminary lattice results show a large spatial b.c. dependence 
for the cc system in the % c \ (J p = 1 + ) channel even near T c , which seems consistent with Ref.[^]. 

6. Summary and conclusions 

We have investigated J and t] c above T c with anisotropic quenched lattice QCD to clarify 
whether the cc systems above T c are compact quasi-bound states or scattering states. We have 
adopted 0(a) -improved Wilson quark action with renormalized anisotropy a s /a t = 4. Anisotropic 
lattice is technically important for the measurement of temporal correlators at high temperatures. 
We have use j8 = 6.10 on 16 3 x (14 - 26) lattices, which correspond to T = (1.11 - 2.07)7;. 
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Table 3: The energy of the cc system in the J channel (J p = 1") in PBC and APBC at J3 = 6.10 and 
p = 0.2fm at each temperature. The statistical errors are smaller than O.OlGeV. We list also uncorrelated 
X 2 /N D¥ and AE =E (APBC) -£(PBC). 

temperature fit range £(PBC) [% 2 /N DF ] £(APBC) [X 2 /N DF ] AE 

l.WTc 7-11 3.05GeV [0.14] 3.09GeV [0.61] 0.04GeV 

1.32r c 8-11 2.95GeV [0.34] 2.98GeV [0.33)] 0.03GeV 

1.6ir c 6-9 2.94GeV [0.10] 2.98GeV [0.22] 0.04GeV 

2.077; 5-7 2.9 lGeV [0.03] 2.93GeV [0.04] 0.02GeV 



Table 4: The energy of the cc system in the r\ c channel (J p = 0~) in PBC and APBC at =6.10 and 
p = 0.2fm at each temperature. The statistical errors are smaller than O.OlGeV. We list also uncorrelated 
X 2 /N DF and AE = £(APBC) -£(PBC). 

temperature fit range £(PBC) [# 2 /#df] £(APBC) [X 2 /N dp ] AE 

l.lir c 7-11 3.03GeV [0.04] 3.02GeV [0.17] -O.OlGeV 

1.32r c 7-11 2.99GeV [0.78] 2.98GeV [0.82] -O.OlGeV 

1.6ir, 6-9 3.00GeV[0.31] 2.97GeV [0.38] -0.03GeV 

2.077; 5-7 3.01GeV [0.03] 3.00GeV [0.07] -O.OlGeV 



To conclude, we have found almost no spatial boundary-condition dependence of the energy 
of the low-lying cc system both in J/*¥ and r\ c channels even on the finite- volume lattice. These 
results indicate that Jf¥ and r\ c survive as compact cc quasi-bound states for T = (1.11 — 2.07)T C . 
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